arXiv:1507.03818vl [math.AP] 14Jul2015 


Generalized gradient flow structure 
of internal energy driven phase fleld systems 


Elena Bonetti 

Dipartimento di Matematica, Universita degli Studi di Pavia 
Via Ferrata, 1, 1-27100 Pavia, Italy 
E-mail: elena.bonetti@unipv.it 

Elisabetta Rocca 

Weierstrass Institute for Applied Analysis and Stochastics 
Mohrenstr. 39, D-10117 Berlin, Germany 

and 

Dipartimento di Matematica, Universita degli Studi di Milano 
Via Saldini 50, 1-20133 Milano, Italy 
E-mail: elisabetta.rocca@wias-berlin.de 
and elisabetta.rocca@unimi.it 

July 15, 2015 

Abstract 

In this paper we introduce a general abstract formulation of a variational 
thermomechanical model, by means of a unified derivation via a generalization 
of the principle of virtual powers for all the variables of the system, including 
the thermal one. In particular, choosing as thermal variable the entropy of the 
system, and as driving functional the internal energy, we get a gradient flow 
structure (in a suitable abstract setting) for the whole nonlinear PDE system. 

We prove a global in time existence of (weak) solutions result for the Cauchy 
problem associated to the abstract PDE system as well as uniqueness in case of 
suitable smoothness assumptions on the functionals. 

Key words: gradient flow, phase field systems, existence of weak solutions, unique¬ 
ness. 

AMS (MOS) subject classification: 74N25, 82B26, 35A01, 35A02. 


1 Introduction 

In this paper we introduce a general derivation of thermo-mechanical phase transition 
models by use of a generalization of the principle of virtual powers, in which micro- 
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forces and thermal forces are included. It is known that a recent held of research, 
in the framework of phase transitions, has concerned models with some micro-forces 
(see, e.g., the approaches by Fremond m and by Gurtin ini)- The main idea is that 
the equations governing the evolution of phase transition phenomena may be derived 
by a variational principle, i.e. the principle of virtual powers, in which micro-forces, 
responsible for phase transitions (i.e. for changes in the microstructure level of the 
materials), are included. As a consequence, the resulting PDE system provides an 
intrinsic variational structure, at least concerning equations for displacements and 
internal quantities, as phase or order parameters. Many authors have dealt with this 
kind of approach. We mention, among the others we quote some contributions as |T5| , 
do], and dl]- 

On the other hand, as far as thermal properties are concerned, in the recent 
years several efforts have been spent to investigate models in which an entropy balance 
(or imbalance) equation was introduced in place of the more classical “heat equation”. 
We recall, e.g., the contribution by 12], i, and [5]. In particular, let us mention 
that the last paper shows a derivation of the equation on the entropy by convex 
analysis tools and the application of a Legendre transformation for the free energy. 
It is interesting to observe that in this framework, also thermal memory is formally 
justihed from the point of view of the derivation of the model. In a different direction 
Podio-Guidugli, in relation to a theory proposed by Green and Naghdi, introduced 
the possibility of including thermal displacements and forces in the whole balance 
of the principle of virtual powers, so that the entropy equation may be recovered, 
as well as the momentum equation, as a “balance of forces”, forcing the system on 
the base of some “reluctance to order”. Indeed, starting from the consideration that 
some virtual power principle may be used to deduce all balance and imbalance laws of 
thermomechanics, he suggested to use it also for the derivation of thermal evolution, 
through the notion of thermal displacement. As a consequence, he derives an equation 
for the entropy of the system, which is combined with momentum balance. This 
approach turns out to be consistent with thermodynamical principles. See, among 
the others, m and the papers by Green and Naghdi na, 1121, and references therein. 
Finally, we can quote the recent contribution [16], where a gradient structure of systems 
in thermoplasticity is introduced by means of a free entropy functional instead of the 
internal energy, which is the driving functional in the present contribution. 

Indeed, in this paper, we aim to combine the previous approaches and provide 
a general abstract formulation of a variational thermomechanical model which can be 
applied to recover different phase transitions and phase separation phenomena, also 
accounting for mechanical or thermal effects. Hence, we introduce an unihed approach 
which formally justihes the evolution of the thermal variable (represented here by the 
entropy of the system), the phase parameters, and (possibly) the displacements. Ac¬ 
tually, in the following, we are dealing with two state variables: s, which mainly plays 
the role of the entropy, and a phase parameter x, representing the internal mechanical 
variable. The main advantage of the gradient structure is the possibility of deriving a 
time-incremental minimization procedure, where the internal energy functional is min¬ 
imized with respect to the entropy and the internal variables and so the existence of 
weak solutions for the associated Gauchy problem can be deduced under quite general 
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assumptions on the involved nonlinearities. 

Indeed, the choice of the energy functional and the dissipation potential are 
fairly general. In particular, in the internal energy functional we can include mul¬ 
tivalued operators to ensure some internal constraints. Since the resulting gradient 
flow structure is nonlinear and non smooth, we have to introduce a suitable notion of 
(weak) solution in order to get a global in time existence result. However, the weak 
notion of solution we are introducing is naturally in accordance with the physical 
meaning of the problem under consideration as well as with the classical principles of 
thermodynamics. The proof is performed by means of a combined regularization and 
time discretization procedure. Moreover, uniqueness of solutions is proved under some 
further smoothness assumption on the internal energy functional. 

The paper is organized as follows. In the next Section [2] we derive the model 
and state the main assumptions on the involved physical quantities and functionals. 
The main existence result is stated and proved in Section [3l as well as the uniqueness 
of solutions. 


2 The model and the main assumptions 

Let H c be a bounded and sufficiently regular domain with boundary T := dfl. We 
introduce an Hilbert triplet H c iL c W (with dense and compact injections), where 
V = H = T^(H), and H is identihed as usual with its dual. We introduce the 

notations (■, ■) for the duality pairing between V and V and (■, ■) for the usual scalar 
product both in H and in To simplify the notation, we write H in place of 

L^(H)^, or V in place of when vector-valued functions are considered. For 

every f e V we indicate by / the spatial mean of / over H, i.e. 

where |H| stands for the Lebesgue measure of H. We note as Hq, Vq and Vq' the closed 
subspaces of functions (or functionals) having zero mean value in H, V, and in V, 
respectively. Then, by the Poincare-Wirtinger inequality, 

ll^^llvo := |Vu|^dx^ 

represents a norm on Vq which is equivalent to the norm naturally inherited from V. 
In particular || • ||vo is a Hilbert norm associated to a scalar product ((•, •))vb (dehned 
in fl2.ip j. and thus we can introduce the associated Riesz isomorphism mapping A : 
Vq Vq by setting, for u,v e Vq, 

(2.1) (Au,v} := {{u,v))vo ■= Vu-Vvdx, 

Jn 

so that {Au,u} = ||u||^, for every u e Vq and <(u,H“^(u)) = ||u||y/ for every v e Vq. 
Finally, we can identify Hq with H'q by means of the scalar product of H so to obtain 


4 


Energy driven systems 


the Hilbert triplet Vq Hq cz Vq, where inclusions are continuous and dense. In 
particular, if z e E and f e Vq; h is easy to see that 

(2.2) Vz ■V{A~^v)(ix = {z — z)v(ix= zvdx. 

Jn Jn Jq 

In what follows in this section we introduce our modelling approach and the set of 
PDEs and initial and boundary conditions which we are going to analyze in the next 
sections. 

2.1 The Principle of Virtual Powers 

The model is derived by using a variational principle in mechanics which is known 
as (generalized) principle of virtual power. Indeed, we refer to some generalization 
of the well known mechanical principle as we are including in the involved forces the 
microscopic forces, acting on some “micro-scale”, and also possible “thermal forces”. 
Without entering the details of this argumentation, let us point out that this principle 
is formally based on the fact that velocities are considered in a suitable linear space 
and thus forces are dehned as elements acting on velocities with respect to some 
duality relation between the two spaces. This is done for any (sufficiently smooth) 
subdomain T> Vt. Hence, before proceeding we make precise the virtual velocities 
we are considering. More precisely, let us take the couple of virtual velocities {St,Vt) 
(whose physical meaning may change time to time). In the case when no accelerations 
are included, the principle of virtual powers can be written considering the power of 
internal forces Vint and of external forces Text (depending on V, St, Vt) as follows: 

Vinti'V, St, Vt) + Vexti'^-i Vt) = 0. 

We assume that the power of internal forces is introduced as follows (in H and for any 
virtual velocities h* e Vb and vt^V) 

(2.3) V^nt = ((F, St)) + ((g, Vt)) = (F, St)+ \ Bvtdx+ \ E Vvt dx, 

Jn Jn 

where F, B and E denote interior thermal and mechanical (micro) forces and stresses, 
respectively and the duality relation ((•,•)) is suitably dehned between forces and 
velocities spaces. 

Analogously, the power of external forces is 

Vext=((n,St)) + ((Z,Vt)). 

We let Vt ^ V and St e Vq and we assume there exists Z, z such that 

(2.4) ((TZ, St)) = (TZ, St) and {Z,vt)= [ Z vt + f zvt, 

Jn Jr 

where Z and z stand for the external forces acting in the bulk H and at the boundary 
T, respectively. 
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It is clear that that we are considering a different behavior of the forces on the 
two types of virtual velocities. Indeed, we note that the elements Q and Z are dehned 
as a.e. forces living in the bulk and on the boundary (with suitable summability), 
while we take F and TZ as general as possible to include all the different (and less 
regular) situations we will face. In particular, as it will be clear once we will make 
a precise choice of the actual velocities (cf. fl2.5p and fl2.6p ). of the energy functional 
fl2.1ip . and of the dissipation potential fl2.14p . we aim to write down an equation for 
the thermal variable of conservative type: it will result indeed as a conservation of 
energy, while the equation for the mechanical variable will be on non conservative 
type. This mainly motivates the choice we have made for the power of internal and 
external forces. Other choices are possible (cf., e.g.. Remark 12.8p . but we prefer not 
to move in this direction in the present contribution. 


2.2 The constitutive relations and the PDEs 

The state variables. We are dealing with a physical system governed by the state 
variables (s, x, Vx) whose evolution is ruled by different thermomechanical relations. 
Note that we are distinguishing between different dependence of the energy with re¬ 
spect to the two variables s and x^ we consider, in particular, the gradient Vx but 
not Vs as state variable (cf. fl2.1ip i. This corresponds to the specific choice we have 
done for the forces F, Q and IZ, Z we have made in fl2.3p and 02.41) . 

In order to get the evolution of s, we take the actual velocities as 
where ^ Tg and Vt = 0, in order to get (for 7^ = 0) (F, = 0 for all ^ Vq 

and so we obtain 

(2.5) A{F) = 0 in V^. 

The evolution of x is obtained by integrating by parts in Vint and choosing Z = z = t) 
and hi = 0 as well; 

(2.6) B-divE = t) in W, 
with the no-mass flux through the boundary of R: 

(2.7) E ■ n = 0 on T , 

where we have denoted by n the outward unit normal vector to T. Notice that we 
have chosen here to have 0 ezternal forces TZ and Z only for simplicity of notation. 

Remark 2.1. In the following, we mainly refer to the variable x as a phase or order 
parameter, i.e. related to the micro-structure of the physical system. However, let us 
point out that we could formally include in our procedure the derivation of the (more) 
classical momentum balance equation (letting, e.g., x stand for displacements). In 
this case, the force B has to be equal to 0, due to the principle of rigid motions. 
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The functionals and the main assumptions. We introduce two functionals gov¬ 
erning the evolution and the equilibrium of our (thermo)mechanical system. These 
functionals depend on the state variables and on the dissipative variables, respectively. 
As far as the equilibrium, it is governed by an energy functional, and we choose to 
make use of an internal energy functional (in place of the free energy functional). This 
choice is motivated by the fact that we may interprete s as the entropy of the system 
(see [TB] for a physical justihcation). However, it is well known that, under suitable 
assumptions, the internal energy may be introduced as the Legendre transformed of 
the free energy. 

Before we make precise the choices of the internal energy functional and of the 
dissipation potential, let us introduce a function W, depending on y, as the sum of a 
convex possibly non-smooth part and non-convex but regular function and it satishes 
some smoothness and growth assumptions, in particular, we need: 

Hypothesis 2.2. Assume W{x) = Pix) + 7 (x); where 

(wl) (3 : dom(/5) ^ [0, -l-oo] is convex, proper, and lower semicontinuous, 

(w2) 7 : R ^ M is a function on R. 

Remark 2.3. Notice that particularly meanigful choices of W used in the literature 
of phase transitions (when x assumes the meaning of phase variable) are, for example, 

1. the double well potential W{x) = {x^ — 1)^ 

2. the logarithmic potential W(x) = xlog(x) + (1 “ x) log(l ~ x) ~ x'^ 

3. the double obstacle potential W{x) = -^[o,i](x) ~ where /[ 07 ] denotes the 
indicator function of the interval [ 0 , 1 ] and it is defined as ^[ 0 , 1 ] (a^) = 0 if xe [ 0 , 1 ] 
and /[ 07 ](x) = -l-oo otherwise. 

Moreover, we introduce a function j{6,x) : R x R —> [0, -l-oo] such that 

6 1 -^ j{9, x) is a convex, proper, and lower semicontinuous for every y e R and 

X j (^) x) is a function for every 6* e R, 


and let 




(2.8) 


if {9,x)^HxH 

and j(0,y) e L^{n) 


[+00 

if {9,x)^HxH 

and j(0,y) ^ ^^(H) 

(2.9) 

Jvi9,x) = Jh{9,x) on 

V X H. 



Hence, we can introduce the convex conjugate of Jy as follows Jy(s,x) : V' x H 
[0, -l-oo] is dehned as 

X) = sup «s, e) - Jv{9, x)) , (s, y) e W X H. 

eeV 


( 2 . 10 ) 
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Now, we are in the position of introducing the energy functional e : V' x H x H 
(- 00 , +oo]: 


(2.11) e(5, X, Vx) = X) + (^1 Vxr + W{x)^ dx. 

Let us note here that the hrst term in fl2.11|] contains both the purely caloric part of 
the energy functional (i.e. the one depending only on s as well as the coupling terms 
depending on both s and x) (cf- Subsection l2.3l for possible choices of Jy), while inside 
the integral over we have the parts accounting for the nonlocal interfacial energy 
effects (the | VxP) and the mixing potential W (cf. Remark l2.3l for examples of possible 
choices of functions W). We intentionally choose not to consider interfacial (nonlocal) 
energy effects in the variable s in order to differentiate the roles of the caloric and the 
mechanical parts (s and x, respectively) in our approach. 

Then, we can dehne the subdifferential (with respect to the variable s) dv',vJv(s, x) 
which maps V x H into 2'^ as (cf., e.g., [6]): 

( 2 . 12 ) 

V e dviyJy(s, x) if and only ii v eV, 

(■s, x) e T'(Jy), and Jy(s, x) ^ <s - w, v) + Jy(w, x) V(w, x) ^ V' x H . 

Actually, in what follows we will always work in the space V x H and so we 
will state directly the assumptions we need on the functional Jy dehned in fl2.inp . In 
particular, we need the following assumptions: 

Hypothesis 2.4. We assume that Jy : V xH ^ [0, +oo] is such that: there exist two 
positive constants Ci,C 2 e such that the functional Jy defined in fl2.10p satisfies: 

(ji) X Jy(s, x) is Frechet differentiable in H for every s e V, 

^ ciWpWh + C 2 , for every r] e dv'yJv(s,x) and (s,x) e D(J^), 

H 

where ^ denotes the partial derivative with respect to x (which will be denoted also 
by d^(-) and by (•)^ in the paper). 

Moreover, we assume that 

(J3) s I—» Jy(s, x) is proper, convex and lower semicontinuous from V to [0, +oo], for 
every x^ H, 

so that the subdifferential dv'yJv(s,x) which maps V x H into 2^ according to the 
definition 02.121) turns out to be a maximal monotone operator acting from V' to 2^, 
for every x^ H (cf |I|j. 

Note that, the assumption (J3) follows from assumptions on j and 02.101) and 
that possible examples of functions j complying with our assumptions will be listed 
in the next Subsection 12.31 


(J2) 


SJv(s,x) 


dx 
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Remark 2.5. Observe that the assumptions on the positivity of the maps j and 
Jy could be weakened: we need indeed to have only a lower bound (possibly with 
a negative constant) for them in order to perform the first a-priori estimate fl3.16j) . 
Moreover, let us nothe that the assumption (J2) could be relaxed: we could indeed 
assume Ci, C 2 to be two continuous functions of x bounded on dom(/5). However, we 
put ourselves in this setting to avoid further technicalities for the reader’s convenience. 

Remark 2.6. Note that, under particular assumptions on the function j (for example 
in case dom(j) = Mj, we could also rewrite the functional e in 02.111) as (cf Remar k \2. 71 
for more details and a) 

(2.13) e(s,x,Vx) = £(r(s,x) + ^|Vxr + l^(x)) dx, 

where j* is the conjugate function of j with respect to the variable s, i.e. 

J*(s, x) = sup (sO - j{d, x)) , V(s, x) e ffi X R . 

SgR 


We introduce as dissipative variables the time derivatives St and Xt (see, e.g. 
|llj . for a definition of the pseudo-potential of dissipation a la Moreau) and we include 
dissipation in the model by choosing the following form for the pseudopotential of 
dissipation depending on the dissipative variables St and Xt- Note that we suppose the 
evolution to be rate dependent. The hrst possibility we consider for the dissipation 
functional is 

(2.14) :p:VjxH^R, ip{st,Xt) = ^(st, ^ f \xt\^ dx . 

Note that we have a natural scalar product in Vj ((•,-))vp' dehned as ((st,Sf))yj = 
{st,A-^St) = ^^Ixtpda; = {xtWt) = Ixifn- In the dehntions of e and (p we 

have normalized all the physical constants to 1 for simplicity and without any loss of 
generality. Another possibility consists in letting 

(2.15) <p{st, Xt) = A-ht) + A-^Xt), 

but we prefer not to exploit this case in the present contribution in order to distinguish 
between the roles of the two variables: the thermal variable s {conserved) and the 
mechanical variable x (non conserved). Moreover, a rate-independent model could be 
introduced in place of the rate-dependent one we analyze here by suitably modifying 
the choice of the dissipation functional fl2.14p (cf., e.g., [IZ]). However, the analysis 
we are performing does not apply to this case, which would require ad hoc techniques 
and some suitable notion of weak solution. 

The constitutive relations and the PDEs. Now, according to the dehnition of 
Vint and of e and ip (cf. 02.111) and 02.141) 1. we let the thermal force F be 


( 2 . 16 ) 


F = d,e + d„<fi = dv’,vJv(s,x) + A '(s,). 
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Hence, for the evolution of x we prescribe the following mechanical (micro) forces and 
stresses B and E\ 

(2.17) B = d^e + = d^Jy{s, x) + d^ix) + x'ix) + Xt, E = = Vx • 


From fl2.5H2.7p and the above constitutive relations we deduce the following 
PDF system for the evolution of s and y: 

(2.18) St + Ar] = Oin Vq, f] e dv',vJv{s,x)y a.e. in( 0 ,r), 

(2.19) Xt - Ax +^ + 7 (x) + d^J*(s,x) = 0, ^ e/3(x) a.e. in 12 x (0, T) 

(2.20) Vx ■ n = 0 a.e. on P x (0, T ), 


where we denote by the subdiferential of f3 {(3 = dj3) and by 7 = 7 '. 

Notice that system fl2.18^ - 112.19p can be rewritten in terms of the vector u : = 
{s, x) in a more general framework, as the gradient-flow associated to the functional 

(2.21) «F(n) = «F(s, x) = x) + £ I Vxl' + W^(x)) dx , 

as 

( 2 . 22 ) 




norm 


Afiut) -f — 9 0 in (0,T), 
ou 

where Af is the duality map between y. := Vq x H and Vq x H induced by the 

(2.23) ||m||^ := <A"^(s),s> + r \x? dx, so that A/'(s, x) := (^"^(s), x) • 

Jq 

Remark 2.7. Let us notice that in case dom(j) = R, which is also equivalent to 
assume 

j*{r) 


lim 


= -Foo , 


|r|^+oo |r| 

then we can prove that (cf. 121 ) the functional e defined in fl2.1ip can be rewritten as 

e(s,x, Vx) = (^f{s,x) + ^|VxP + W{x^ dx, 

where j* is the conjugate function of j with respect to the variable s, i.e. 


J*(s, x) = sup (s0 - j(0, x)), V(s, x) e R X R . 

deR 

Moreover, the inclusion fl2.18p can be rewritten as the following gradient flow in Vj: 

St + dy/Jy(s, x) 3 0 , 

where dydy E defined as the subdifferential of Jy in V mapping V' x H into 2^' as 
follows: 

(2.24) 

^ e dv'Jy(s, x) if and only if ^ e V', 

{s,x)^D{J^), and Jy{s,x) ^ is - +Jy{w,x) V(n;, x) e R'x R , 
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where ((•,•))* denotes the scalar product in V. The reader can refer to and to 
Section 2] for the proofs of these results. Finally in this case we have u e dv'yJv{s, x) 
in V iff u E dsj*{s,x) a.e. in where dg denotes here the subdifferential of convex 
analysis with respect to the variable s (cf, e.g., l^). 

2.3 Possible choices of j* 

In this section we show how to derive different types of phase-field models by our 
general system. 

2 2 

The Caginalp model of phase transitions. Choose j*{s,x) = ^— sy + 
Denote hy 9 := dgj* = s — x- Then, the Hyp. 12.41 is obviously satisfied and the PDEs 
fl2.18f[2H^ can be rewritten as 

6'i xt - A6' = 0 , 

Xi - Ax + ^(x) + 7(x) -^30, 

coupled with Neumann homogeneous boundary conditions on 6 and x- This PDE 
system can be easily identified with the “standard” phase field model of Caginalp 
type (cf. [7j), letting 9 be the relative temperature of the system and x the local 
proportion of one of the two phases of the substance undergoing phase transitions. 

The entropy model for phase transitions. Choosing j*{s,x) = j*{s — A(x)) = 
exp(s — A(x)), we have that Hyp. 12.41 is satisfied in case A is a Lipschitz continuous 
function on the domain of f3. Then, defining 9 := dgj* = exp(s — A(x)), we get 
s = log9 + A(x) and the PDEs 02.18412.15]) can be rewritted as 

(log 0 -t-A(x))t-A 0 = 0 , 

Xt - Ax + /3(x) + 7(x) - A'(x)^ 3 0 , 

again with Neumann homogeneous boundary conditions for both 9 and x- This system 
can be easily identified with the “entropy” phase field model introduced in [3] and 
[1]. Here 9 respresents the absolute temperature of the system which is forced to be 
positive, by the presence of the logarithmic nonlinearity in the 0-equation. Let us 
notice that in this case the assumption D{j*) = M is not verified, hence we are not 
entitled to use the function j* instead of the operator Jy in e (cf. Remark l2.6p . so, 
the choice we made here is only formal. For a rigorous analysis of this case the reader 
can refer to [5]. 

The Penrose-Fife model for phase transitions. We choose j*{s, x) = ~ log('S — 
x), for s > X and define 9 := Aj*(s,x) = “ 7 ^- Then, we observe that we can 
formally get the Penrose-Fife mode. Indeed, it results that d^j*{s,x) = 7 ^ = |- 
Thus, we can rewrite fl2.18f[2H9l) as 

(» + x)i +aAI) - 0 , 

Xt - Ax + /?(x) + 7 (x) + ^ ^ 0 > 
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coupled with Neumann homogeneous boundary conditions on — ^ and y. This system 
can be easily identihed with the Penrose-Fife model of phase transitions introduced in 

ng. 

Remark 2.8. Let us notice that in case we choose as pseudopotential of dissipation 
the functional fl2.15p in fl2.16p and fl2.17lh the first equation is the same as fl2.18lh 
while the equation for x results as: 

(2.25) A-\t -Ax + I3{x)+ lix) + x) 9 0 , 
and thus 

(2.26) st + A x)) 3 0 in , 

(2.27) Xt + Aw = 0, we-Ax + /3(x) + 7(x) + ^xf(^’X), 

coupled with Neumann homogeneous boundary conditions for dsj*(s,x), X w. 
In this case the evolution of x L ruled by the well-known fouth order Cahn-Hilliard 
equation modelling phase separation phenomena (cf, e.g., l^). However, as we already 
mentioned, we prefer not to deal with this case here. 


3 Main results 

In this section we state the main results of this paper, the hrst one (Thm. 13.11) concerns 
the existence of global in time solutions for system fl2.18p - fl2.19p coupled with the 
boundary condition fl2.7p and the initial conditions 

(3.1) s(0) = So in D(J^), 

(3.2) ^(0) = xo a.e. inR, 

while the second one (Thm. 13.41) regards uniqueness of solutions under more restrictive 
assumptions on the nonlinearities involved. Let us start with the hrst result. 

Theorem 3.1. Assume Hvpotheses \2. 41 and \2.2\ and take sq e D{Jy), Xo ^ V ndom(^). 
Then, for every T > 0 there exists at least one solution (s, x) to fl2.18p - fl2.20p and 
fl3.1l) - fl3.2p satisfying the regularity properties: 

(3.3) ssH\0,T;V'), 

(3.4) X e H\0, T; H) n L”(0, T; V) n ^^(o, T; H^Q)) . 


Proof. In order to prove Theorem l3.ll we hrst approximate system fl2.18f[2T^ 
with a regularized problem depending on a positive small parameter e and then we 
pass to the limit by (weak-strong) compactness arguments and semicontinuity results 
based on sufficient a-priori estimates - independent of e - we are going to prove on 
the approximating solutions. 
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The approximated problem. Let us fix e: > 0. Then, for every T > 0 and 
(so,e, Xo) ^ {D{Jy)r^H) X [V ndom(^)), we aim to find a solution (sg, Xe) ^ 

H) to the following differential inclusions: 

(3.5) dtSe + A{r]e + ese) = 0 in Ej, rje e dv’,vJv{se,Xe), for a.e. t e (0,T), 

(3.6) dtXe-^Xe +ie + l{Xe) + d^Jy{se,Xe) = ^ ie^l^iXe), a.e. in X (0, T) , 

coupled with the boundary and initial conditions fl2.20p and fl3.1fl3.2p . with so,e in 
place of Sq. In particular, we assume that 

(3.7) So,£ e 7)(J^) n H, so,£ ^ Sq in Ej as e \ 0. 


We first observe that we can recover fl3.5p and and fl3.6p . by approximating the 
energy functional 02.211) as follows: 

®£(s,x) := ^£(s,x) + £ + ^{x)^ ^'^(sw) := <7v(s,x) + | £ Is^dx. 

Actually, note that now is defined in (Eq r\ H) x H. Hence, we can construct its 
subdifferential in the duality between Vq and Vq, and rewrite the equation 03.5p as 

(3.8) dtSe + AQ = 0 in Eq, Ce ^ Sv',v^s{se, Xe) for a.e. t e (0, T ). 

Now, our aim is to prove the existence of solutions of 03.5p - 03.6p . 03.ip - 03.2p with so,e 
istead of Sq, and 02.2Up by a time-discrete approximation, as follows (cf. also |2T] for 
a similar procedure). Here we drop the index e in order to simplify the notation. Let 
us fix a time step t = T/N, N e N and introduce a uniform partition 

Pr := {to = 0, ti = T, . .., tn = nr, ..., tN = T} 

of the interval (0,r). Then, we need to find a discrete approximation s"' ~ s{tn), 
X" ^ x(tn) by solving the implicit Euler scheme (cf. also fl2.22p L 

/Tj^ — 

(3.9) Afl - - - j+C” = 0, n = l, N-U^ ■.= uo, 

where C"" ^ and we have defined = (s"',x"'), Uq = (so,£)Xo)- Using the 

functional space, we have already introduced to define the operator JV, T-L = (Eq n 
H) X H, we notice that fl3.9p is the Euler equation for the variational problem 


j find [/” e 77 minimizing 

j FeiT,U^-\U):=j^\\U-U^-% + <l>,{U), UeH. 

It is not difficult to see that this minimization problem is solvable due to the lower- 
semicontinuity and coercivity properties of $£ (cf., e.g., |211|22] and references therein 
for a similar variational approach to find a discrete solution). 

Then, we can construct the piecewise constant interpolants Urit) '■= if 
t E [[n — l)r, nr]. In particular, we get that and we recover the solution U : = 





















Bonetti - Rocca 


13 


(s, x) of fl3.5H3.6p as the limit of [ 7 ^ as r \ O.This can be done, by nsing snitable a 
priori estimates (independent of r and then passing to the limit by compactness and 
semicontinnity argnments. We do not enter the details of the proof, as it is very similar 
to the estimates and passage to the limit procednre we are going to detail in the next 
sections to pass to the limit as £ \ 0. Note that in this case some technicalities are 
avoided dne to the more regnlar setting for the variable (recall the strict positivity 
of e). Thns, we can easily prove the following existence resnlt. 

Theorem 3.2. Under the same assumptions ofTheorem lS.ll letting e > 0 behxed and 
03.71) holds, then there exists a solution to 03.5^ - 03.61) with Se(0) = so.e and Xe(0) = xo, 
with the following regularity 

(3.11) Ss e H\0, T; Rj) n L®(0, T; R), 

(3.12) X, e H\0, T; H) n L®(0, T; V) n T; 


A priori estimates (uniform in e). Let us consider the system 03.5H3.6I) . where 
for the sake of simplifying notation we neglect the index e for solutions and involved 
functions.We now perform the a priori estimates independent of e, so, we use here the 
same symbol c for positive constants, possibily different from line to line, depending 
on the problem data, but independent of £. 

In order to perform the hrst estimate we need to prove here a small variant of 
the chain rule formula stated, e.g., in [9l Prop. 4.2]. 

Proposition 3.3. Let G : V' x H ^ [0, -l-oo] be a map such that 

V G{u, v) is Frechet differentiable for every u e V , 

u G{u, v) is a proper convex lower semicontinuous mapping for every v e H, 

and let u e H\0,T;V') n L\0,T;Vo), v e H\0,T; H) n L^{0,T;V), and S{t) e 
dv'yG{u{t),v{t)) for a.e. t e (0,r), where the subdifferential dy'y is dehned as in 
fl2.r2p . Then the function g = G{u{-),v{-)) is absolutely continuous in [0,T] and 
g'{t) = (u\t),5{ty) + {v\t),dvG{u{t),v{t))) for a.e. fe(0,T). 

Proof. Here we follow the lines of the proof of [9l Prop. 4.2]. Let w e 
ppi, 00 ( 0 , 2 ^) non-negative function with compact support in (0,T). We choose 

h > 0 such that supp(tc) c [h,T — h]. For a.e. t e (0,T), by dehnition of sub¬ 
differentials we can infer that 

(u{t) — u(t — h), S{t — h)} + {v{t) — v{t — h), dyG{u(t — h),v{t — h)) 

^ hif) ~ d(f — h) < (u{t) — u{t — h), S{t)} + {v{t) — v{t — h), d^G{u{t), v{t)). 

Indeed, observe that (5, dyG) belongs to dG, i.e. to the sub-differential of the function 
G : V' X H [0, -foo] dehned w.r.t. the variable {u,v). Observe that we can extend 
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w outside of (0,T) with the 0 value. Hence, multiplying by w{t), integrating with 
respect to t, and letting h \ 0, we obtain 


h 


— (u{t) — u{t — h), — h)')w(t)dt = — {u(t + h) — u(t), S(t)}w(t + h)dt 


h 


■T-h 


\ Wit), 

Jo 


6it)}wit)dt, 


h 


T idit) — git — h))wit)dt = — 5it)iwit) — wit + h))dt 


h 


—> — r 5it)w'it)dt, 

Jo 

— r (uit) — uit — h), Sit))wit)dt ^ r Wit), <^it)Wit)dt. 
h Jh Jo 


Moreover, 


h 


h 


r (u(t) — — h), dyGiuit — h),vit — h))wit)dt 

Jh 

1 

" T ivit + h) — vit),dyGiuit),vit))wit + h)dt ^ iv'it),djjGiuit),vit))wit)dt, 

h Jo Jo 

I (u(t) — — h), —> I iv'it),d^Giuit),vit))wit)dt. 

Jh Jo 


Therefore, we conclude that 

rT rT 


— r git)w'it)dt = r i(Wit),dit)') + iv'it),d„Giuit),vit))))wit)dt 

Jo Jo 

for all non-negative Lipschitz continuous test functions w with compact support. Since 
both the positive and the negative part of a Lipschitz continuous function are Lipschitz 
continuous, we obtain the assertion. I 

First a priori estimate. We test fl3.5p by A~^st getting 


(3-13) — ||st||y^/, 

and in addition, using the dehnition of Jy and a variant of the chain rule formula 
stated in Proposition 13.31 with the choices G = Jy, u = s, u = x, we get 

(3.14) {Ag, A-^st) = (st, g) = ^Wi^it), xiW - ixt, ^xWi^, x)) ■ 

Testing fl3.6p by Xt, we get 

(3-15) ^ f ^W) + iXt,SxJvi^,x)) = tJ ■ 

Adding the resulting equations and integrating over (0,t), t e (0,r), and using the 
dehnition of we obtain 

(3.16) (llstllvg'+ IImIIh) dT + Ts(s(t),x(t))+ £hE(x(^)) + l|Vx(f)P ^ c. 
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where here c depends in particular on the initial data. Adding to both sides in fl3.16p 

\\x{t)\\H = llxolli + 2 f (x(r),xt(r)) dr, 

Jo 

and using Holder and Young inequalities together with Hyp. I2.2l and [2^ and a standard 
Gronwall lemma, we obtain 

(3-17) l|■S^llL2(0,T;V'g') + ^ll■5|lL«>(0,T;/^■) ^ 

(3-18) \\x\\H^{0,T-,H)nL'^{0,T-,V) ^ C, 

(3-19) ||/^(x)||l®(0,T;L1(O)) ^ C- 

Second a priori estimate. We proceed by a comparison in fl3.8p . Due to fl3.17P i . we 
have that is bounded in L^(0,T; Vq), and thus 

(3-20) ||C£||L2(o,r;Vo) ^■ 

Hence, using fl3.17l) o. we get 

(3-21) \\Ve\\L^{0,T;H) ^ C. 

Indeed, due to the dehnition of 'hg we can infer that Cs = Ve + where we have 

used the notation | s^dx and the fact that has sa domain the whole 

real line and thus its subdifferential in the duality V, V corresponds to the standard 
sub differential of the convex analysis (cf. [6]). 

Third a priori estimate. Using now Hyp. 12.41 together with fl3.2ip . we get 
(3-22) X)\\l'2{0,T;H) ^ Ci||he||r2(0,T;ir) + C 2 < c. 

Moreover, by comparison in fl3.6p and by standard monotonicity and regularity results, 
we get 

(3.23) ||'C||l2(o,t;H) + \\x\\L^{o,T-,m{n)) ^ c. 


Passage to the limit as e \ 0. Now, we aim to pass to the limit in fl3.5p - fl3.6p 
as e \ 0, revovering dually a solution to fl2.18112.lOp . By virtue of fl3.17h3.23p and by 
compactness results we get (at least for some subsequences of £ \ 0): 


(3.24) 

* 

Se^S 

in H\0,T-X), 

(3.25) 

Xe^X 

in H\0, T- H) n L®(0, T; V) n L2(0, T; i7^(D)), 

(3.26) 

Ce-C 

in L\0,T-Vo), 

(3.27) 

Ve^V 

in L2(0,T; H) , 

(3.28) 


,Xs)^z in L2(0,T;i7) , 

(3.29) 


■0 in L®(0,r;R) . 
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Notice that, by the dehnition of H/g and by fl3.26p . fl3.27p . and 03.281) . we immediately 
deduce that C = r] a.e.. Moreover, by strong compactness, from fl3.25p we can also 
deduce (at least) (cf. [23]) 


(3.30) 


mC\[0,T];H)nL\0,T;V) 


We aim to identify C e dyiyJy{s,x) (see 03.261) 1. By dehnition of sub differential this 
corresponds to prove that 


(3.31) 


<n-s,C>dr<r (Jy(n,x) - Jy(s,x)) dr Vn e Eg' andx^H. 


Note that, if we test the equation 03.Sp by A and integrate in time, by weak lower 

semicontinuity of norm we have (for e dv',v^e{se, Xe)) 


(3.32) 


I 


■T nT 

hmsup<S£, C£>dr ^ <s,Odr. 
0 e\0 Jo 


Hence, by employing 03.26p . 03.27p . and 03.29p . and using the fact that (e belongs to 
Xe),we get, for all v e Eg and e H, 


(3.33) 


<'i'-s,C>dr^ <n,C>dr - limsup (se,QdT 
Jo Jo £\0 Jo 

T + liminf (-(se, Ce)) dr 

£\0 Jg 

< liminf f ((n, Q - Ce)) dr 

£\0 Jg 


< limsup (d'£(n,Xe) - d'£(s£,Xe)) dr 
e\0 Jo 

= — lim inf 

e\0 

^ — lim inf 

£\0 


Jg i~'^viv,Xe) + Jvise,Xs) - ^\\v\\h + dr 

i-JU'f^^Xe) + Jv{Se,Xe)) dr. 

Using Hyp. 12.41 and 03.3Up . we can deduce 

(3.34) 

Hence, let us observe that 

(3.35) J(^(5„ X.) - X) = JviSe, Xe) " J^^e, X) + x) - J^(s, x), 
where 


lim J*(n,X£) = -/f(t,x) ■ 
£-\0 


and 


h-.= JUSe,Xe)-JviSe,X): 
h ■= JU^e,x) - Jyi^^x)- 
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Hence, by fl3.22p and fl3.3Up . we have 


(3.36) r |/i|dr^r \\d^Jy{se,Xe)\\H\\Xe- x\\h<^t as £ \ 0 . 

Jo Jo 

The second integral 

r / 2 dr=r - Jy(s,x)) dr 

Jo Jo 

is treated by using the fact that Jy is lower semicontinuous in with respect to weak 
convergence as it is convex (for y hxed) and so 


■I 


lim inf | J 2 dr ^ 0 

e\0 


Hence, coming back to 03.331) . and using 02.4p . we have, for all v s Vq and e H, 


liminf {-Jy{v, Xe) + Xe)) dr 

0 e\0 

T rT 


e\0 


hniinf(-J^(n,ye)) dr - linr inf (J^(s£, y^) - Jy(s, y) + Jy(s, y)) dr 




lim inf R 

Q £\.0 


0 £\0 

T rT 

lim inf I 2 

Q £\,0 


f •/v(s,x)dr 

Jo 


< f {JyiP^x) - Jvi^^x)) dr, 

Jo 


and this concludes the proof of 03.3ip . Finally, using the previous convergences, where 
we have now identihed ( = r] e dv',vJy{s,x) H), we can pass to the limit 

in the approximated system 03.5H3.6I) as well as in the corresponding boundary and 
initial conditions for e \ 0. This concludes the proof of Theorem 13.II I 


We conclude now with the last result of our paper concerning uniqueness of 
solutions for problem fl2.18fl2.20p . fl3.1fl3.2p . 

Theorem 3.4. Assume Hvpotheses \2.4 and \2.2\ and take sq e D{Jy), Xo^V ndom(^) 
and suppose moreover that 

(3.37) the maps y t—> dyyJy^s, y) and y t—>. d^Jy{s, y) are Lipschitz continuous 

from H to Vq and from H to H, respectively, for every s e Vh Then the solution 
(s, y) of li2.1t^fT2ni) . lid. 103.21) is uniquely determined and the following continuous 
dependence estimate holds true: 

(3.38) 

||(si - S 2 )(t)||y, + ||(yi - X 2 ){t)\\H + |V(yi - X 2 )\\HdT (||(si - S2)(0)||^, 

+ ll(xi - X2)(0)||^) . 
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Proof. Let us take the difference of the two equations 02.181) and the two 
relations 02.19P corresponding to two different solutions {si,Xi), i = 1,2 and test 
them by — S 2 ) and (xi — X 2 ), respectively. Integrating over (0, t), for t e [0, T], 

using Hyp. 12.41 and Hyp. 12.21 we get 



Using then the Lipschitz continuity of 7 and assumption 03.37p . together with Gronwall 


lemma, we obtain exactly 03.381) . 
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